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Abstract 

Using the method of shape invariant potentials, a number of exact solutions of one di- 
mensional effective mass Schrodinger equation are obtained. The solutions with equi-spaced 
spectrum are discussed in detail. 

1 Introduction 

Since the introduction of the factorization method by Shrodinger to solve the hydrogen atom 
problem algebraically, there has been a considerable effort to generalize this method and to find 
exactly (analytically) solvable potentials for the Schrodinger equation. One of the pioneering 
works in this line was done by Infeld and Hull 2 . They obtained a wide class of exactly 
solvable potentials using the factorization method. After the introduction of supersymmetric 
quantum mechanics :SJ> the concept of a shape invariant potential (SIP) was introduced by 
Gendenshtein 4. in 1983. The Schrodinger equation with a SIP has exact solution and its 
energy eigenvalues can be obtained algebraically. The method of shape invariant potentials is 
also applied to some other equations like Fokker - Plank equation as well [Hj . Although SIPs are 
widely studied and many families of them are obtained fSj , but the problem of the classification 
of SIPs has not been completely solved yet. 

In this paper we apply the method of SIPs to effective mass Schrodinger equation and obtain 
its exact solutions in a specific ansatz. Effective mass Schrodinger equation was introduced 
by BenDaniel and Duke [2] in 1966 to explain the behavior of electrons in the junctions of 
semiconductors. It also have some applications in heterostructures, graded alloys, and quantum 
wells 0, U, EH- 

Recently, effective mass Schrodinger equation is studied by some authors and some of its 

exact solutions are obtained |12j . |13| . Here we show that the method of shape invariant 

potentials applied to this problem leads to a wide class of exact solutions. 
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The paper is organized as follows. In section 2 we review the method of SIPs and its 
generalized form for effective mass Schrodinger equation. In sections 3 and 4 we study two main 
classes of generalized shape invariant potentials(GSIPs) corresponding to special choices for the 
parameter dependence of the potential. Such parameter dependence have been considered to 
obtain SIPs for ordinary Schrodinger equation in reference [5] and references therein. Due to the 
importance of quantum systems with equi-spaced energy levels, we investigate them separately 
in each case. Section 5 is devoted to summary and concluding remarks. 



2 Generalized Shape Invariant Potentials (GSIPs) 

In this section we review the method of SIPs for effective mass Schrodinger equation. Consider 
an eigenvalue problem given by 

Htp n {x) = E n ip n (x) , (1) 

where H is given by 

H = -4-tA^4- + V(x). (2) 

dx2m(x)dx v 1 K ' 

Here m(x) is the position depended effective mass and V{x) is the potential. 

Suppose that the minimum eigenvalue of H is Eq. Let us define a new Hamiltonian Hi := H— 
Eq. The Hamiltonian H\ can be written in the following form 

in which V\(x) = V(x) — Eq. The ground state energy of H\ is zero, i. e. the spectrum of H\ 
is non negative. Consequently the Hamiltonian H\ can be assumed to be 

#i = A* A , (4) 

where A is a first order differential operator of the form 

A = U(x)4- + W(x) . (5) 
dx 

Inserting A from Eq. (JSJ) in Eq. (j3J and comparing the result with Eq. © one gets 

1 



, , -U 2 (x), (6) 
V x {x) = W 2 (x) - [U{x)W{x)]' . (7) 

Here, a "' indicates the derivative with respect to x. Eq. Q is a generalization of the well known 
Riccati equation. We call the solution of this equation [3], W(x), the superpotential. 
Next we define the Hamiltonian H2 as follows 

Then the potential V%{x) reads 

V 2 (x) = W 2 (x) + [U(x)W(x)}' - 2U'(x)W(x) - U(x)U"(x) , 

W 2 (x) - [U(x)W(x)}' + 2U(x)W'{x) - U(x)U"{x) . (9) 



The Hamiltonian Hi is called the super symmetric partner of H\ It can be easily shown that 
Hi and H2 has the same spectrum except for the ground state of H\. Suppose that ipi(x) is an 
eigenfunction of H% with eigenvalue E\ . Then in view of Eq. © we have 



A^Aipi(x) = E x il)x{x) 



(10) 



Using Eq. (jHJ) it is easily seen that Aipi (x) is an eigenfunction of H2 with eigenvalue E\ , provided 
that ^i(x) is not the ground state wave function. 

One can repeat the above procedure for H2 and get the supersymmetric partner of H2. Then 
one arrives at a hierarchy of Hamiltonians which their spectrum are essentially the same. Now 
suppose that the potentials V\{x) and V^(x) depend on some parameters which we show them 
in a compact form by a. The shape invariance condition is 



V 2 (x,a 1 ) = Vi(x,a 2 ) + -R(ai) , 



(11) 



where a>2 = F(ai) is a function of a\ and R{a\) is independent of x. If these conditions are 
fulfilled, the spectrum of H\ can be found algebraically jSJ: 



E n = Y,R(a n ) 



(12) 



In what follows, we use this method to find the exact solutions of effective mass Schrddinger 
equation considering the super potentials in the ansatz: 



W(x, a) = ag(x) + f(x) + -h(x) . 



(13) 



Inserting W(x,a) from Eq. (|T3"|) in Eqs. Q and @, the shape invariance condition (fTTj) gives 
{al-a\)g 2 {x) + 



h 2 (x) + 2(a 2 - a x )g{x)f{x) + 2 (— - -) h(x)f{x) 



U(x)(a 2 - a x )g{x) + U{x) (— - — ) h(x) 

\a 2 aij 

+U{x)U"{x) + R{a l ) = . 



1 1 



2U(x) 



a\g{x) + —h(x) + f(x) 
ai 



Before starting to study Eq. (|14j) . we make a definition for the future convenience: 

f ' x dx' 



Y x) 



U(x') 



(14) 



(15) 



There is also a comment in order on the square integr ability of the wave functions. According 
to Eq. Q the ground state of Hi is given by Atpo(x) = 0. This, together with definition (JSJ) 
leads to 

W(x' 



ipo(x) ~ exp 
Then the square integrability of V'o ( x ) leads to 



U(x') 



-dx' 



ipQ(x)dx 



exp 



W(x') 
U(x') 



dx 



< 00 . 



(16) 



(17) 



This condition puts some restrictions on acceptable W(x) and U(x). In fact the following 
condition must be fulfilled 

l im fE^L x ' = +00 . (18) 

If has definite signs for x — ► ±00, then there is a simple condition which guarantees 

Eq. (HBJ rrj. Let us define 

— := lim — ±-L . 19 

V U J ± x^±oo JJ{x) V ' 



Then, to fulfill Eq. ifTBJl . it is sufficient that 



,'W\ [W\ . , 

sign — = -sign — . (20) 



17/ + V U J _ 

It can be shown that this condition also ensures the square integrability of excited states wave 
functions [""]• 

3 GSIPs of type W(x,a) = f(x) + a 

In this section we consider GSIPs whose superpotentials are of the form 

W{x,a) = f(x) + a . (21) 

This means that in Eq. 1)13(1 we have put g(x) = 1 and h(x) = 0. Then the shape invariance 
condition gives 

U(x) [U'(x) - 2/0*0]' " («2 - ai) [U'(x) - 2f(x)) + (a 2 2 - a\) + fl(oi) = . (22) 

In this equation C/(x) and f(x) are independent of parameter a±. Therefore there are only 
two ways for Eq. (|2*2*)) to be consistent. The first way is that U'(x) — 2f(x) = const. =: 2uq. 
Then the superpotential is given by W(x,a) = ^U'(x) + a — uq. and Eq. ((2*2*]) implies that 
R(a±) = {a\ — a|) — 2uo(ai — 02). Absorbing uo in a, one gets i?(ai) = [a\ — a|). Therefore 
in this case 02 is an arbitrary function of ai provided R(a\) = a\ — a| > to ensure that the 
energy eigenvalues are non negative. 

The second way is that a\ — a 2 and (a| — af) + i?(ai) be independent of parameter, i. e. 

R(ai) = a(2 ai -a) + R , (23) 
a 2 = ai - a , (24) 

where a and i?o are some real constants. Using these conditions and defining u{x) := U'{x) — 
2f{x) we can write Eq. ((22(1 in the following form 

U(x)u'(x) + au(x) +Ro = . (25) 

Before trying to solve this equation let us have a look at the spectrum of H\. The energy 
eigenvalues of H\ are given by 

n 

E n = Y^ R(ai) = n(R + 2a x a - no?) . (26) 

i=l 



This equation shows that for a = the energy levels are equi-spaced. Due to the importance 
of systems with equi-spaced energy levels we consider two distinct cases corresponding to a = 
and separately. We call solutions with equi-spaced energy levels oscillator like solutions. 



3.1 Oscillator like Solutions 
For a = Eq. ([25]) gives 

u(x) = -RqY(x) , (27) 
where Y(x) is defined by Eq. (|15[). Therefore the superpotential is 

W(x, a) = X - [U'(x) + R Y(x)] + a , (28) 

and the potential Vi(x,ai) is given by 

Vtix, a) = - [(RoY(x) + 2a) 2 - 2R ] + V (x) , (29) 



in which 



V (x) := -\u'\x) - \u{x)U'\x) . (30) 



The energy eigenvalues are E n = uRq ,n = 0, 1,2, ••• . Since the energy eigenvalues are non 
negative we should have Rq > 0. because the energy eigenvalues are non negative. 
Now we can obtain the zero energy eigenstate using Eq. (|16|) 

^o(x) ~ |[/(x)|~ 1/2 exp (^-^Y 2 (x) -aY(x) \ . (31) 

It can be easily verified that the ground state wave function is normalizable i. e. Eq. (|17|) is 
fulfilled. The results of this subsection is summarized in tabic ^ 

3.1.1 Examples 

Example 1: If U (x) = ^= = const, the system is a harmonic oscillator; the most simple system 
with equi-spaced energy levels. 

Example 2: For U(x) = 5- we have y(x) = x 2 . The effective mass and potential are given by 

m{ps) = 2x 2 , (32) 

Example 3: For C/(x) = eo& 1 h - we have y(x) = sinhx and m(x) = ^cosh 2 x. This form of 
position dependent mass is considered by some authors in graded alloys jlUj . The potential is 
given by 

(34) 



Vx(x,a) = - 



,„ , , 9 sinh x — 2 „ 

(Ro sinhx + 2a) 1 -\ 2R C 

cosh x 



3.2 Exponential Solutions 

If a ^ 0, by redefinition of the function f(x) and the parameter a we can put Rq = 0. The 
redefinition is given by the following equations 

a — ► a = a + ^ , (35) 
2a 

f{x) - /(*) = /(*)-§*. (36) 
Under this redefinition, W(x,a) does not change but R(a±) changes as 

R(ai) — ► -R(ai) = a(2ai — a) , (37) 
and Eq. (|25|) can be rewritten as 

U(x)u (x) + au{x) = . (38) 
We conclude that for a ^ one can always put Ro = 0. The solution of Eq. (|3*S|) is 

u(x) = uq exp (— aY(x)) . (39) 
We call these solutions, exponential solutions. The superpotential W(x, a) is therefore 

W(x, a) = - [U'(x) - n exp (-aY(x))] + a , (40) 
and the potential is given by 

V^x, a) = ^ [(n e- ay(a;) - 2a) 2 - 2au e~ Qy(:c) ] + Vb(x) , (41) 
in which Vq{x) is defined in Eq. ()3U|) . The energy eigenvalues are 

S n = an{2ai - an) . (42) 
The ground state which is given by Eq. Q16j) takes the following form 

Mx) ~ ^(x)!" 1 / 2 exp _ ai y( x )^ . (43) 

The above results are summarized in the tabled 

Examples 

Example 1: U(x) = Uq = const. With this choice for U(x) we arrive at the Morse potential for 

ordinary Schrodinger equation which is a well known SIP [3]. 

Example 2: U(x) = —fx. In this case, choosing no = 1 in Eq. Q39j) . one gets 

W(x,a) = -- (a + 2x 2 ) +a , (44) 

1 a 2 
Vi (x, a) = -x 4 - [a/2 + a)x 2 - — + a 2 , (45) 
4 lb 

2a l 1 / x 2 \ 
i>o{x) ~ |x|~"2 exp ( -— j . (46) 



4 GSIPs of Type W(x, a) = ag(x) + f(x) 

This ansatz is obtained putting h(x) = in Eq. (|13|) , Applying the shape invariance condition 
we arrive at 

[a\ - a\)g 2 (x) + (a 2 - ai )g(x) [2f(x) - U'(x)] 

-(01 + a 2 )U(x)g\x) - U(x) [2f(x) - U'(x)]' + R(a x ) = . (47) 

Assuming that a\ — a 2 = a, one gets 

a 2 g 2 {x) - 2af(x)g(x) + a [U(x)g(x)}' - 2U{x)f(x) + U(x)U"(x) 

-2ai {ag 2 (x) + U{x)g'{x)) + R( ai ) = . (48) 

Since U(x), /(x), g(x), and a are parameter independent, we conclude from the above equation 
that 

a 2 g 2 (x) - 2af(x)g(x) + a [U(x)g(x)]' - 2U(x)f'(x) + U(x)U"(x) + C x = , (49) 
-2ai (ac/ 2 (x) + U(x)g'(x)) + i?(ai) - Ci = . (50) 

In the same way Eq. ()5UI) implies that 

ag 2 (x) + U{x)g' (x) + C 2 = , (51) 

fl( 0l ) = Ci - 2aiC 2 . (52) 

Therefore energy eigenvalues are 

E n = n 2 aC 2 + n[d- {2a x + a)C 2 ] ■ (53) 

Using Eq. (j5T|l in Eq. one finds 

a5(x)n(x) + U(x)u'(x) + C\ - aC 2 = , (54) 

where u(x) = U'(x) — 2f(x). Eq. (|5^|) indicates that for aC 2 = the energy eigenvalues are 
equi-spaced. In the following we discuss three different families of solutions corresponding to 
aC 2 = 0, aC 2 > 0, and aC 2 < 0. 

4.1 Oscillator like solutions (aC 2 = 0) 

These are solutions with equi-spaced spectrum (see Eq. 1)53(1 ) . Three different cases can be 
realized as follows, 

Case 1: a = C 2 = 0. In this case according to Eq. (|51|). g(x) should be a constant function 
of x. We have discussed this case in section 3. 

Case 2: a = and C 2 ^ 0. In this case Eqs. (|51|) and (|54l) are easily solved and we get 
g(x) = —C 2 Y{x) and u(x) = —C\Y{x). Therefore, 

W(x, a) = hrf{x) + \(Ci- 2aC 2 )Y{x) , (55) 
E n = n((Ci - 2aiC 2 ) . (56) 



One can redefine the parameter oasa-»a=y - aC 2 - The results are summarized in the 
table El 

Case 3: a ^ and C 2 = 0. In this case from Eq. (|5T|) one finds 

= «kr) ■ (57) 

Then substituting this solution in Eq. (|54jl we have 

u(x) = -^y(x)--^L, (58) 
2 y (x) 

where C3 is a constant. In fact C3 is not so important and can be absorbed in the parameter a. 
See table H 

4.2 Trigonometric Solutions (aC 2 > 0) 



Putting k := \JC 2 a it is straightforward to solve Eqs. (joTj) and (|54j) to get 

k 

g (x) = tan [JfeYYsc)] , (59) 

a 

and 

_ Ci - aC 2 _ ^ ^ _ 

Then the superpotential is obtained: 

W(x, a) = (-— + Cl 7^ ) tan(fcy) + — sec(A;y) . (61) 
a 2k 2 



Redefining the parameter a as 



/eft Ci k , . 



we have a := d\ — a~ 2 = —k. Using this redefinition one can write the general form of superpo- 
tential as 

W(x, a) = -a tan(«y) + b sec(aY) + -U'{x) , (63) 
where b = — #. The energy eigenvalues are given by 

E n = na(na — 2a\) . (64) 

Finally the ground state is given by 

ip (x) ~ \U(x)\^-\sec(aY)\^\tan(aY) +8ec(aY)\-« . (65) 



It is remarkable that for U(x) = const, the above solution reduces to the Scarf I poten- 
tial which is a SIP for ordinary Schrodinger equatioin [Sj. The results of this subsection is 
summarized in table [5] 



4.3 Hyperbolic Solutions (aC 2 < 0) 



Considering the definition k := \J — C201 one can solve Eqs. (|5Tj) and (f51|) to get 



and 



k 

g(x) = - tanh [kY(x)\ , (66) 
a 

= _ Cl ~ a ° 2 tanh fifey(x)] - C 3 sech \kY(x)] . (67) 



Then the superpotential can be found as 



W(x, a) = (— + Cl 1^ ) tanh(£y) + ^sech(fcy) . (68) 
a 2fc 2 

Now we redefine the parameter a as 

~ ka Ci + k 2 
a^a = — + — — — . (69 
a Ik 

After this redefinition of a we have a := d\— 6,2 = k. Therefore the general form of superpotential 
can be given as 

W(x, a) = atanh(aY) + teech(aY) + ^U'{x) , (70) 
where b = The energy eigenvalues are given by 

E n = na(na — 2a\) , (71) 

and the ground state is 

ip (x) ~ |f/(x)|"2 1 |sech(ay)|^|tanh(ay) +sech(ay)|~- . (72) 

For U(x) = const, the above solution reduces to the Scarf II potential which is a SIP for 
ordinary Schrodinger equatioin 0. The results of this subsection are summarized in table 

5 Concluding Remarks 

In this paper we studied exact solutions of effective mass Schrodinger equation using the mathod 
of shape invariant potentials. We considered an ansatz in which the superpotential takes the 
form 

W(x,a) = ag(x) + f(x) + ^-. (73) 

a 

Assuming h(x) be vanishing, the exact solutions for the Schrodinger equation are studied and a 
wide class of solutions with equi-spaced spectrum are obtained. Following the general arguments 
of the paper, it can be verified that if h(x) is non-vanishing then it should be simply a constant. 
In this case the shape invariance condition requires m{x) to be constant. Thus if h(x) 7^ 0, one 
obtains the SIPs for the ordinary Scrodinger equation. The main results of the this study are 
summarized in tables 1-6. 
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Superpotential 


W(x, a) = \ \U'(x) + RqY(x)] + a 


Potential 


Vi(x,a) = \ [(R Y(x) + 2a) 2 - 2R ] + V (x) 


Energy eigenvalues 


E n = tiRq 


Ground state 


ifo(x) = N\U{x)\~ 1 ' 2 exp (-^Y 2 (x) - aY(x)) 



Table 1: GSIPs of type W(x, a) = f(x) + a; Oscillator like solutions 



Superpotential 


W(x, a) = \ [U'(x) - u exp (-aY(x))] + a 


Potential 


Vi(x,a) = \ [(u e- aY W - 2a) 2 - 2au Q e~ aY ^] + V (x) 


Energy eigenvalues 


E n = an(2a\ — an) 


Ground state 


ifoix) ~ \U{x)\ V2 e ^(x) aiYC*) 



Table 2: GSIPs of type W(x,a) = f(x) + a; Exponential solutions 



Superpotential 


W(x,a) = \U'(x) + aY{x) 


Potential 


Vi{x, a) = a 2 Y 2 (x) - a + V (x) 


Spectrum 


E n = 2a\n 


Ground state 


Mx) ~ \u{x)\~ 1 ' 2 exp(^y 2 (x)) 



Table 3: GSIPs of type W(x,a) = ag(x) + f(x); Oscillator like solutions (I) 



Superpotential 


W(x,a) = lU'(x) + ^Y(x) + ^ 


Potential 


V^x, a) = § Y 2 (x) + (g + §) vki + ^ - 2) + F ° (X) 


Spectrum 


£"„ = Cm 


Ground state 


^o(x) ~ |C/(x)|-V2|^ (a ,)|-f exp(-f y 2 (*)) 



Table 4: GSIPs of type W(a;,a) = ag(x) + /(x); Oscillator like solutions (II) 



Superpotential 


W(x,a) = -atan(ay) +bsec(aY) + \U'{x) 


Potential 


Vi(x,a) = (atan(«y) — 6sec(ay)) 2 

— asec(aF)(6tan(ay) — asec(ay)) + Vq{x) 


Spectrum 


E n = na(na — 2a\) 


Ground state 


V'o(^) ~ \U(x) \~\ sec(aY)\~\ tan(ay) + sec(aY)\~& 



Table 5: GSIPs of type W(x,a) = ag(x) + f(x); Trigonometric solutions 



Superpotential 


W(x,a) = atanh(ay) + bsech(aY) + \U'(x) 


Potential 


Vi(x,a) = (atanh(ay) - fcech(ay)) 2 
+asech(ay)(6tanh(ay) — asech(oy)) + Vq(x) 


Spectrum 


E n = na(na — 2a\) 


Ground state 


ip (x) ~ |f/(a;)|^ 1 |sech(ay)|~o|tanh(ay) +sech(ay)|~« 



Table 6: GSIPs of type W(x,a) = ag(x) + f(x); Hyperbolic solutions 



